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1. Introduction 

There are two kinds of billiards in (n + 1) dimensional unit cube D n+1 , 
(n > 1), differing by accepted reflection law, are considered in the classical 
mechanics and dynamical systems theory. It is the usual billiards with the 
reflection law - the angle of incidence is equal to angle of reflection, and the 
bouncing billiards, when a particle, moving whithin cube rectilinearly and 
uniformly in a given direction 9 = (9 1 , . . . , 6 n+l ), if reaching any edge of 
cube is bounced vertically to the opposite one and then goes parallelly to its 
previous path. On the other hand, there are known two different kinds of one 
dimensional periodical motion: the rotation and libration [1-4]. In the case 
of bouncing billiards each of the coordinates of the orbit's point performs 
the rotation, whereas in the case of usual billiards each of them performs the 
libration; in either cases the frequencies are equal to components of vector 9. 
The orbits of both motions, if considering on the (n + 1) dimensional torus 
7 i or ( ri + 1 ) anc } cu be D n+1 respectively, are the continuous lines which either 
are closed (periodical motion) or, under assumption that 9 is nonresonance 
vector, according to H. Weyl theorem ([5], see also [2, 3, 6]), are distributed 
uniformly on the whole surface of Tor^ n+1 ^ or in cube D n+1 (conditionally 
periodic motion). For instance, the Kepler problem on motion of planets [2] 
can be considered in terms of bouncing billiards, while for the problems on 
small vibrations of the pendulum [1, 2], the Lissazhu figures [1, 2], and the 
vibrations of atom [1] it is natural to consider the orbits are given in the 
cube, in cartesian coordinates, and not on the torus. 

We obtain the terms Xk G D n of the (Z-th, 1 < I < n+ 1) billiard sequence 
from the consecuitive points Zk of intersections of an usual billiard's orbit with 
two opposite edges of D n+ , if erase from every vector Zk the (Z-th) component 
which is equal to or 1. The main result of present paper establishes that 
under some simple assumption on the reduced frequency 9 l+1 /9 l , the (z-th, 
1 < % < n) by-coordinate billiard sequence (x^ > ) k % 1 can be splitted into 
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countable set of countable pairwise disjoint subsequences in such a way that 
the restriction of the return map 7Tj : x$ — > x^ +l on each of them appears 
the rotation transformation of the form 

a + b{a + pk} (A; = 1,2,...) (1) 

of the circle {x : \x — a\ = \b\}. Such kind of splitting into finite num- 
ber of subsequences with the simultaneous exact estimating the density of 
arising exceptional sets is also considered. We note that the return map 
7T = (7Ti, . . . , 7r n ) is similar to Poincare map, which inherits [4] the basic 
features of the whole orbit. 

It is known that the libration occurs in the neighborhood of the elliptic 
points in phase space [1-4, 7]. Hence, in spite of our Theorems 1-4 relate to 
billiards in cube, we note that through use of canonical transformations [1- 
3] , they are also applicable to the integrable and weak perturbated Hamilton 
systems. By this reason, the corresponding theoretical-mechanical reformu- 
lations - the Corollaries 1-3, are also provided. 

Note in the end, that due to the results below, the coexistence in weak 
chaotic Hamilton systems of two (undistinguished in small bins of phase 
space) different kinds of motion, mentioned by I. Prigogine (for the case of 
weak stable systems, [7], Ch. 2), is reduced to known coexistence of the dis- 
crete and continuous forms. 

2. Main results 

This work contains improved formulations of our results on billiard boundary 
sequences, announced before in the note [8]; the corresponding theoretical- 
mechanical reformulations are also given. In order to emphasize the underly- 
ing dynamical context, in the latter statements the terminology "condition- 
ally periodic libration in cube" is used. In the "act ion- angle" coordinates 
[1-4, 9-11], if we consider the "angle" variables as the cartesian coordinates, 
the orbits of such motion are exactly the usual billiard's orbits. 

As in [8], we define the billiard sequence of usual billiard in cube D n+1 
as follows. Let for a pregiven index 1 < I < n + 1 the Zk be all the points of 
an orbit for each of which = 0, 1 (it is supposed that the orbits do not 
intersect the ribs of cube D n+1 and the index k > 1 increases with time). 
The sequence Bg = (x k )^ =l of those points from D n , for which x^ = 
whenever 1 < s < I — 1 and x^ = z^ +V> whenever / < s < n we call the 
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(7-th) billiard (boundary) sequence. We note that the Poincare section of 
bouncing billiard's orbit consists of the points of the form (see [4]), 

({«i + fck}, ...,{«„ + Pnk}) (A; = 1,2,.. .) (2) 

({.} is the fractional part of number). It is not difficult to show, for usual 
billiards their billiard sequences Bg are of the form 

(\\a 1 + p 1 k\\,...,\\a n + P n k\\) (A: = 1,2,...) (3) 

where the quantity 

| \x\ | =2 min({a;}, 1 — {x}) 

is so-called "the distance to the nearest integer" of point x G R 1 , that is 
used in the theory of Diophantine approximations [6, 12, 13]. In contrast to 
sequences of fractional parts from Eq. (2), which are studied quite well (see 
details in [6]), the properties of sequences (3) are almost unknown (unless the 
results from [8], [14] and several interesting but isolated number-theoretical 
relations, can be found in [6], Ch. 2). 

Our results relate to structure and properties of the by-coordinate se- 
quences Bg i = (x^)^ =1 and their return map 7Tj : x$ — > x^ +l under two 
assumptions, which are mentioned in all further formulations: 

l/2<Wi<l where = |0 (i+1) /0 (,) | (1 < i < n) , (4) 

and the (reduced) frequency oji is irrational number. The restriction (4) 
(which is equivalent to such one, accepted in [8]) on the billiards under con- 
sideration is arised from the method applied. This condition appears to be 
very natural and is quite enough (due to some property of symmetry of the 
discrepancy) in order to research the discrepancies of the general sequences 
of fractional parts ([14]). 

The next Theorem 1 establishes the property of splittability of billiard 
by-coordinate sequences Bg i into countable number of subsequences of frac- 
tional parts of form (1). Theorem 2 estimates the rate of the split process in 
dependence on arithmetical characteristics of the frequency a;*. The metrical 
Theorem 3 and the Theorem 4 estimate this rate. Also, we formulate the 
theoretical-mechanical analogies of these theorems, the Corollaries 1-3, are 
based on representation of by-coordinate return map 7Tj in a form of (count- 
able or finite) set of the discrete rotation transformations of circles. 
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Theorem 1 If the condition (4) is satisfied, then there exists a partition of 
the natural series N into countable number of countable and pairwise disjoint 
subsets N ijS , 

oo 

N = Y,N hS , N hSi nN hS2 = (si^s 2 ), (5) 

s=l 

such that for some numbers a ijS , b itS , a i>s , (3 i>s and every s > 1 the equalities 

x 2 s (fc) = a i,s + Ksi a i,s + Pi,* k } (A; = 1,2,...) (6) 

hold. Here the function v i s is numbered the set N i>s in increasing order of 
its terms and each /5 ijS is irrational of the form /5j )S = T iiS (u;j) where r,^ s is a 
fractional-linear transformation with entire coefficients. 

What is more, it can be also shown, that for each of the coefficients /3i iS from 
Eq. (6), one of the following 5 simple relationships on equivalence, 

/3i, a ~Wi, (3i, s ~u}i/2, f3 itS /2~Ui, 2f3, hS ~u h 2(3 i>s ~ 1 - Ui 

holds. We recall (see [12, 13]), that for large enough values of index k the k- 
th coefficients of expansion into regular continued fraction of two equivalent 
numbers are coincided. 

In the formulations of Corollaries 1-3, we use the terms "rate of split- 
ting" and "boundary motion". The first one means the density of the corre- 
sponding exceptional set, while the latter means the motion of terms of the 
by-coordinate billiard sequences = determined by return map 

tT = (ttx, . . . , 7r n ) and ir^x^) = Tr^ji) (1 < i < n; k = 1, 2, . . .) . 

(7) 

Indeed, if we introduce a discrete time t! = {1,2, . . .} in such a way that 
the fall of moving particle on two opposite edges of D n+1 occures at its 
consecuitive moments (k and k + 1), the Eq. (7) can be understood as the 
equation of motion of the point x£ in respect of t'. Furthermore, through 
the function v i>s defined in Theorem 1, we can introduce the splitting of the 
time series t' into countable number of discrete independent times t' i s {= N), 

k e t' i s <=> u i>8 (k) e N itS . 
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(i) 

Then, according to Theorem 1, the motion of point x v , fc s in respect of 
time t' i s = {k : k > 1} appears a uniform discrete rotation of the circle 
{x : \x — <ii tS \ = \bi jS \} with the frequency /5j )S . The function i/ ijS can also be 
explicitly expressed by some characteristics (the functions A and r below) of 
the frequency a;,. 

Corollary 1 Lei us have a conditionally periodic libration in cube D n+1 and 
for some 1 < % < n the condition (4) is satisfied. Then the boundary motion 
of i-th coordinate of point Xk G Bg can be splitted into a countable set of 
discrete uniform rotations. 

In order to formulate the next two theorems we need some definitions. For 
numbers x G (0,1/2) we consider the expansion into continued fraction of 
the form (see 



x 



e o| , e l| e k 



\a\ \a 2 \a>k+i 



+ t-^- + ■■■ (8) 



where eo = 1 and for k > 1 the coefficients = a,k(x) and = €k(x) are 
defined as follows 

/ [l/^- 1 ] {l/^ 1 } < 1/2 / 1 {l/^- 1 } <l/2 

afc ~\ [1/V>H + ! {l/^ 1 }>l/2 ' efc ~j-l {l/^ fe ~ 1 }>l/2. 

Here [.] and {.} are the entire and fractional parts of number respectively, 
ip°(x) = x, and for k > 1 ip k is fc-th iterate of the function ip, 

1>(x) = min ({l/x}, 1 - {l/x}) (=^||l/x||). (9) 

We use also the binary coefficients pk- Pk = whenever e& = — 1 and pk — 1 
whenever = 1 and rewrite the expansion (8) in the form 

x= [a 1 ,a 2 ,...;p 1 ,p 2 ,...\ ■ (10) 

This modified expansion appears to be more appropriate tool for consider- 
ing problems (as well as for problems concerning the discrepancy, see [14]), 
than the regular continued fractions. There is the next interrelation between 
the coefficients ak(x) and Pk(x) and the coefficients bk(x) of expansion into 
regular continued fraction of number x (see [8]): 

a-k = 6 ff (fc) + 2- p k - p k -!, p k = <=> b aik )+i = 1 (11) 
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where the quantity a (= a x ) is defined as follows: cr(l) = 1, a(k + 1) = 
a{k)+2 — pk (k > 1). We define also the functions A (= A x ) and r (= t x ), are 
mentioned in our formulations: A(l) = 1, X(k + 1) = X(k) + r(k) + 1 where 
t(1) = 1—pi and r(fc + l) = r{k) whenever a\(k+i) i s odd, r(k + l) = l — r(k) 
whenever a\(k+i) is even. For a set E C AT its density is defined as 

dens(E) = \im^p~ 1 \{k E E :1 <k <p}\) (here |.| is power of the set) . 

The Theorems 2-4 express the rate of the split process, established in 
Theorem 1 (and which, we have assumed, is measured by the density of 
exceptional set), by the number-theoretical characteristics of the reduced 
frequency cjj. 

Theorem 2 // the condition (4) is satisfied, and the function \i (= /ij) is 

M = 2p + 2 jr Pm+1 (p=l,2,...) (12) 

fc=i, T (fc)=i 

where A = A w ' , t = t u >. , and oo^ = 1 — Ui , then for every p > 1 there is a 
finite partition of natural series into pairwise disjoint subsets, 

W=E N hS + E hP (13) 

s=l 

suc/i i/ioi /or eac/i 1 < s < // and all k > 1 the equations (6) hold and for 
the density of exceptional set E ijP we have 

A(p+1)-1 A(p+1) 

den S (^, P ) = II ^ s (^') = II M«0> • • • ; • • •] • ( 14 ) 

s=0 s=l 

By virtue of the relationships (11), the value of dens(E itP ) in (14) can also be 
expressed directly through the coefficients bk of the frequency In such a 
way, we have all the quantities A, p and dens(Ei iP ) in explicit form. Hence, it 
is not difficult to obtain from Theorem 2 many interesting particular relations 
giving the dependence of density of exceptional set E itP on the arithmetical 
properties of the coefficients and pk or bk- 

The next statement is a reformulation of Theorem 2. 
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Corollary 2 Let us have a conditionally periodic motion in cube D n+1 and 
for some 1 < i < n the condition (4) holds. Then for every p > 1 the 
boundary motion of i-th coordinate of the point x k G B^ can be splitted 
into a finite set of fa{p) discrete uniform rotations and the split rate is equal 
to right-hand expression in Eq. (14)- 

Theorem 3 // the condition (4) is satisfied, then for a.e. (in sense of 
Lebesgue measure) numbers uji e (1/2,1) and enough large p > 1 the se- 
quence i can be splitted into fa, 3p < fa < 4p subsequences of fractional 
parts of the form (6) in such a way, that for the density of exceptional set 
E itP we have 

dens(E ijP ) = (1 + o(l)) exp (-7 p) (p — > 00) 
where < 7 < 00 is an absolute constant. 

By use of the explicit form of invariant measure of the ergodic shift ip, 

V>([ai,a 2 , . . . ;pi,p 2 , • • •]) = [02,03, • • • 5 P2, P3, ■ ■ •] , (15) 

defined by Eq. (9) (see [8]), one can obtain the exact numerical value of this 
constant (see [8]): 

3 L* 2 (^I)-L2 2 (±^) 
7 = 4 ; r V 4 ; = 1.524 . . . 

where Li 2 is the Euler dilogarithm ([15]). It is interesting to note that such 
a constant, is equal to |7, is arised in a metrical theorem on estimating the 
discrepancy (see [14]). 

The point (b) of the next theorem follows from the previous Theorem 3 
and is given again just for completness. 

Theorem 4 If the condition (4) is satisfied, then the sequence can be 
splitted into fa 2p < fi < 4p subsequences of the form (6) in such a way, that 
for the density of the exceptional set E itP from Eq. (13), the next statements 
are true: 

(a) there exists a constant C > such that for every Ui e (1/2, 1) we have 
dens(E i>p ) = <3(exp(-Cp)) (p -> 00) ; 
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(b) for almost each Ui G (1/2, 1) we have 

dens{E itP ) = (1 + o(l)) exp(~7p) (p -> oo) ; 

fc^ i/ £ne cjj G (1/2, 1) zs a Liouville number, then we have 

dens(E ijP ) = o(exp(-C(p)p)) (p -> oo) 

where C(p) — > +oo as p — > oo; m dependens on u>i the quantity C(p) can 
tend to +oo as quickly as we like. 

The next corollary is the reformulation of Theorem 4. The distance between 
two motions, as it is mentioned above, should be understood in sense of 
density of the exceptional sets E ijP from Eq. (13). 

Corollary 3 If the condition (4) is satisfied, then the next 3 statements are 
true: (a) the (boundary motion of) conditionally periodic libration is ever 
approximated by a finite set of discrete rotations at least with the exponential 
rate; (b) the (boundary motion of) almost each conditionally periodic libra- 
tion is approximated by a finite set of discrete rotations with the exponential 
rate; (c) there exist the conditionally periodic librations (the boundary mo- 
tions of) which are as close to a finite set of discrete rotations, as we like. 
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